STABILITY OF CONTACT DISCONTINUITY FOR STEADY EULER 

SYSTEM IN INFINITE DUCT 
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Abstract. In this paper, we prove stability of contact discontinuities for full Euler system. 
We fix a fiat duct Ao of infinite length in R 2 with width Wo , and consider two uniform subsonic 
flow = {uf,0 

> Pl> Pi) with different horizontal velocity in Ao divided by a flat contact 
^ discontinuity Y c d- And, we slightly perturb the boundary of A/o so that the width of the 

perturbed duct converges to Wo + w for \ui\ < S at x = oo for some S > 0. Then, we prove 



< 

' of Ao and 5 are sufficiently small, then there exists unique asymptotic state U^r with a flat 



C3 



that if the asymptotic state at left far field is given by f/, , and if the perturbation of boundary 



1. Introduction 



contact discontinuity r* d at right far field(a; = oo) and unique weak solution U of the Euler 
system so that U consists of two subsonic flow with a contact discontinuity in between, and 
that U converges to Uf 1 and at x = — oo and x — oo respectively. For that purpose, we 
establish piecewise C estimate across a contact discontinuity of a weak solution to Euler system 
(f-) ' depending on the perturbation of <9A/o and 8. 

> ' 

CO 
00 

in 

^ \ Let p, u and p be density, velocity and pressure of flow. Then steady inviscid compressible 

flow is governed by the steady Euler system 

div(pu) = 

div(pu u + pi) = (I : identity matrix) (1-1) 
div(puB) = 
with the Bernoulli's invariant B given by 

B= l -\u? + T ^- (1.2) 
2 1 1 ( 7 -l)p y > 

for ideal poly tropic gas with an adiabatic exponent 7 > 1. The quantity c, given by 

r ™ 

is called the sound speed. The flow type is classified by the Mach number M = — . If M > 1 
then the flow is called supersonic, if M < 1 then it is called subsonic. If M = 1 then the flow is 
called sonic. If the flow is supersonic then the system (jl.ip is hyperbolic, if the flow is subsonic 
then the system (jl.ip becomes a hyperbolic-elliptic mixed system. 

Due to the nonlinearity of the system (11. lj) . one expects that a solution of (jl.ip may contain 
discontinuities such as shocks or contact discontinuities even if a boundary condition is given 
by a smooth function. Such discontinuities can be described through a weak formulation of 
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(jl.ip . A shock and a contact discontinuity are characterized by a normal velocity. While the 
normal velocity is nonzero on a shock, the normal velocity on a contact discontinuity completely 
vanishes. More details are given in Section [2J Because of the difference, one needs different 
schemes to study a shock and a contact discontinuity. For the case of a shock problem, one 
can identify a shock as a graph by using nonzero normal velocity on the shock. Owing to this 
advantage, the stability or instability of various shock phenomena have been investigated in 
many works (see [2], [3], [4], [5], [6], [8], [TT], [13], [15] and references therein). For the case of 
a contact discontinuity, on the other hand, we need a different strategy due to the zero normal 
velocity on the contact discontinuity. For that reason, this subject has been studied in restricted 
regimes yet ([5], [7]). In [5], S. Chen proved the stability of steady Mach reflection configuration 
in a bounded region in M? provided that an appropriate constant pressure is fixed on a cut-off 
boundary. In the Mach reflection configuration, two reflected shocks are separated by a contact 
discontinuity. In [7J, S. Chen and B. Fang proved the conditional stability of a reflection and 
a refraction of shocks occurred when an incident shock hits the interface, which is a contact 
discontinuity, of two different media. But still, the structural stability of a contact discontinuity 
under a general perturbation in unbounded domain is unknown. The main difficulty in study of 
a contact discontinuity is that the states on both sides of a contact discontinuity are unknown, so 
one needs to solve a free boundary problem with both sides of a free boundary to be determined. 
And, this is the main difference from a shock problem of Euler system. 

The study of contact discontinuity is essential to understand a Mach reflection, which is one 
of important but difficult subject. When a vertical incident shock hits an inclined ramp, if the 
incident shock is relatively strong or the angle of the ramp is relatively small, then the incident 
shock is reflected at a point away from the boundary of the ramp, and two reflected shocks 
are formed at the reflection point with a contact discontinuity in between. This phenomenon 
is called Mach reflection, named after Ernst Mach. Also, the Mach reflection for steady Euler 
system can be considered through shock polar analysis(see [5j). It is conjectured that the steady 
Mach reflection in M 2 is structurally stable. In order to prove this conjecture, one needs to prove 
stability of a contact discontinuity along with two reflected shocks. In this paper, we prove 
structural stability of a contact discontinuity in an infinite duct where the flow on both sides of 
the contact discontinuity is subsonic. This is related to the case where two reflected shocks are 
transonic shocks in steady Mach reflection. 

We fix a flat duct Afo of infinite length in M 2 with width of Wq, and consider two uniform 
subsonic flow in A/"o divided by a flat contact discontinuity T c d- Then we perturb the boundary 
of Mo with a small C 1,a function so that the width of the perturbed duct converges to Wq + u 
for \u\ < 5 at x = oo for some small constant 5 > 0. Then we want to show that there exists 
two layers of subsonic flow divided by a contact discontinuity in the perturbed nozzle, and that 
the new contact discontinuity is a small perturbation of r 

It is a new feature that we allow for a perturbed contact discontinuity to converge to different 
asymptotic states at x = ±oo. Since the right asymptotic width Wq + w of a perturbed nozzle 
is not necessarily same as the width Wq at x = — oo, we expect for the asymptotic pressure p^ 
at x = oo to be different from the asymptotic pressure p-oo at x = — oo, and this yields two 
different asymptotic states for a perturbed contact discontinuity at x = ±oo. Details are given 
in Section [2j 

Another interesting aspect is application of the result from [12J to this work. Like a shock 
problem, a contact discontinuity problem is a free boundary problem. But unlike a shock 
problem, the states on both sides of a contact discontinuity T is unknown. Moreover since 
the normal velocity on T is zero, it is not clear how to locate a position of T. So we use 
the Euler-Lagrange transformation to reformulate the contact discontinuity problem as a fixed 
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boundary problem. Then the new difficulty is to find a weak solution of a first order nonlinear 
elliptic-hyperbolic mixed system so that the weak solution satisfies the Rankine-Hugoniot jump 
condition on a fixed boundary. If we can show that the weak solution is piecewise C 1 in two 
subregions divided by the flattened contact discontinuity T, then a simple integration by parts 
shows that the weak solution indeed satisfies the Rankine-Hugoniot jump condition on T. From 
this point of view, we employ the result of [12J to achieve piecewise C 1 regularity of weak 
solution for a nonlinear elliptic- hyperbolic mixed system. According to Theorem 1.1 of [12] . if a 
discontinuity boundary, which is a contact discontinuity in our case, is in C 1 '^ for < f3 < 1, then 
the regularity of a corresponding weak solution to a uniformly elliptic equation is weaker than 
piecewise C 1 ' 13 regularity. But, this may cause a difficulty in the iteration procedure which is 
the way to solve our main problem. Fortunately, the Euler-Lagrange transformation transforms 
a contact discontinuity to a flat boundary which is smooth. Hence, there is no deterioration of 
regularity of weak solutions to elliptic equations in the iteration, and this is an advantage of the 
Euler-Lagrange transformation. 

In Section [2j we compute asymptotic states at far field in a perturbed duct, and use the 
Euler-Lagrange transformation to reformulate our problem as a fixed boundary problem. Then 
we state our main theorems. In Section [31 we establish piecewise C 1 regularity of weak solutions 
to Euler system, and use this estimate to prove the main theorems. The main difficulty would be 
uniform L°° estimate of weak solutions to a class of uniformly elliptic equations in unbounded 
domain especially because we have two different asymptotic states at far field x = — oo and 
x = oo. 



2. Problems and main theorems 

2.1. Asymptotic states at far field. Let p, u and p be the density, velocity and the pressure 

of flow respectively. Then steady flow of compressible polytropic gas is governed by the Euler 
system 

div(pu) = (2.1) 

div(pu <g> u + pl) = (2.2) 

div(puB) = (2.3) 

with the Bernoulli's invariant 

B = -\u\ 2 + , 1P . (2.4) 

2 1 1 ( 7 -l)p y > 

for an adiabatic exponent 7 > 1. We note that if (p, u,p) is in C 1 , and satisfies (|2"3 |) -(|23 |) . then 
it also satisfies the transport equations 

u • V(-^) = and u ■ VB = 0, (2.5) 

and this means that the entropy (4^) and the Bernoulli's invariant B are preserved along each 
streamline in C 1 flow. 

We consider flow in M?. Let u\ be the horizontal component of u, and let 112 be the vertical 
component of u. For an open and connected set Q C M 2 , if U = (u±,U2,p, p) satisfies 



(2.6) 



/ pu-D£= / {pu k u + pe k ) D£= puB •£>£ = () 
Jn Jn Jn 

for any £ G Cq°(0) and k = 1, 2, then U is called a weak solution to the Euler system (|2. l|) - f)2.3[) 
in O. Suppose that is divided into two subsets by a non self- intersecting C 1 curve, and 
that U is C 1 in and C° in Q^. Then one can easily check by integration by parts that U is 
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a weak solution of the Euler system if and only if U satisfies (|2.ip - (12.3p pointwisely in and 
the Rankine- Hug oniot jump conditions (abbreviated as R-H conditions hereafter) 

[pu • n] r = (2.7) 
[p(u ■ n)u + pn] r = (2.8) 
pu ■ n[B] r = (2.9) 

for a unit normal n on V where [F]r is defined by [-F(x)]r := F(x)\-^r — F(x)\-^r for x £ T. By 
(|2.7p . the condition (|2.8p can be rewritten as 

p(u • n)[u • r] r = 0, • n) 2 + p] r = (2.10) 

where r denotes a unit tangential on T. For p > in O, the first condition in (12.10P implies 
either it • n = on T or [u ■ r]r = 0. Suppose that 17 is discontinuous on T. If u ■ n ^ 
and [tt • r]r = hold on T, then T is called a shock. If u • n = and [u ■ r]r 7^ 0, then T is 
called a contact discontinuity. From this, we get the R-H conditions corresponding to a contact 
discontinuity as follows: 

un = on T, [p] r = 0. (2.11) 

In R 2 , we consider a flat duct A/"o = {(x,y) £ R 2 : x £ R, — 1 < y < 1} of infinite length 
and two layers of uniform flow in Mo divided by the line y = with satisfying the following 
properties: 

(i) The velocity and density of top and bottom layers are given by positive constants 
(i^~,0),p^~ and (uj , 0) , pj~ respectively with u+ ^ u^; 

(ii) The pressure of both top and bottom layers is given by a positive constant pi; 

(iii) The top and bottom layers are subsonic flow. In other words, there hold 



11 1 u I i 

' " 1 ' " 1 - yiPi/Pj and c r = yjiPi/Pi- ( 2 - 12 ) 



< 1 and — — < 1 for c[ 



From this, we define a piecewise constant vector U\ by 

TTl \ f (u>f,0, Pl ,p+) =:U+ fory>0 

Ui(x,y) :=■> l _ 1 _ (2.13) 

{{Ui ,u,PhPi ) =■ U t for y < 0. 

Then, TJ\ is a weak solution of the Euler system with a contact discontinuity on the line y = 0. 
Let 77 be a smooth function satisfying 

{0 for x ^ — 1 

, < i]'(x) < 10 for all x£i (2.14) 
1 for x > 1 

Fix two constants uj± with < 1, and let h± £ C 1,a (R) be two functions satisfying 

\\h- + 1 + w_??||^i (]R ) + \\h + - 1 - w+r/Hfl-iQR) < a, (2.15) 

||h_ + l||cfi.«»(R) < cr, || /t + - 1 1| ca,<*(R) < cr, (2.16) 

lim \\h- + 1 + w_77|| c n,a (R u_ jB = D lim ll^+ ~ 1 ~ w +^llc 1 >«(K\hR, J RJ) = (2.17) 
for small cr > to be determined later. To satisfy (|2.16p and (I2.17p . we assume that 

|w±| < cr. (2.18) 

For such functions h±, let us set 

M := {(x,y) £ R 2 : h-(x) < y < h + (x),x £ R}. 
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N is a duct perturbed from Afo by the functions h±. Particularly, the width of iV at x = oo is 
changed to 2 + {u+ + w_) from 2. The goal is to prove that there exists a weak solution of the 



y = i- - -- 



Af 

y = -l- 



Euler system in Af with a contact discontinuity, and that the contact discontinuity is a small 
perturbation of y = for uj± and /i-t satisfying (|2.15p - (|2.18p where the parameter a is chosen 
sufficiently small. This is to achieve structural stability of contact discontinuity of steady Euler 
system even when the width of the duct is changed at far field. On the boundaries y = h±(x) 
of Af, we describe the slip boundary condition 

u-n bd = (2.19) 

where is the inward unit normal of N on y = h±(x). 

For a C 1 function h*(x) satisfying h-(x) < h*(x) < h + (x) for all x in R, suppose that 



U{x,y) = (u,v,p,p)(x,y) 



(u+, v + ,p + , p + )(x, y) loxy>h*(x) 
(u',v~,p-,p-)(x,y) for y < h*(x) 



is a weak solution of the Euler system in N with a contact discontinuity on y = h*(x) and u > 
in Af, and that U is C 1 in Af \ {y = h*(x)}, and satisfies the slip boundary condition (|2.19p 
on dAf. Also, suppose that U converges to XJ\ in L°° at x = — oo. Then the total mass flux on 
x = xq in each of top and bottom layers is preserved as well as the total mass flux on x = xq in 
./V is preserved for all xq in R. In other words, U satisfies 

h*(x) ph+{x) 

(p~u~)(x,y) dy = pffuff =: m Q , / (p + u + )(x, y) dy = =: mj, (2.20) 

/ (/w)(a;, y) = pfuf + pfuf =: m . (2.21) 

(l2~20]l and (l2~2T]l can be easily checked by using ([21]), (l2TTT]l and (l2TT9j) . Moreover, by ([23]) and 
positivity of u, U also satisfies 



for 5 = ^ and S defined by (fL2l) . 

By (|2.17p . the boundary of iV gets flat at far field. From this, we expect that the flow in N at 
far field becomes two layers of uniform flow divided by a flat contact discontinuity. So we first 
compute the asymptotic state at far field in Af consisting of two layers of uniform flow with a 
flat contact discontinuity in between. At x = — oo, we fix U\ in (|2.13p as the asymptotic state. 
It remains to compute the asymptotic state at x = oo corresponding to U\. 

Lemma 2.1. Fix 7 > 1, and let U[ be as in ()2. 13j) . Set oj := w+ +uj- for oj± from (j2. 17j) . Then 
there is a constant ojq > depending on Ui and 7 such that for any uj G [— ujq,ujq\, there exist 
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unique constants ,u~ , p+ , p~ , p r and G (— 1 — w_, 1 + so that 

Ur{ X ,v) = if t ^ pt \ =:U T ; {° ry>UJ * (2.23) 
[ (u r , 0,Pr,p r ) =: U r for y < a;* 

is a weak solution of (jl.ip , and satisfies (|2.1ip and (|2.20p - (|2.22p with a contact discontinuity 
y = . Furthermore, we have 

dp r 

-q^ > Po > for all to G [-wo,wo] 
where (3q is a constant depending only on t/j,7 and ojq. 

Proof. One can easily check that U r in (I2.23P is a weak solution of (ll.ip . and satisfies (12. lip . 
Set U r (x,y) = (u r ,0,p r , p r )(x,y). By (|2.22|) . we can write (u r ,p r ) as follows: 



(u r ,p r )(x,y) 



for y > lu* 



2(^- 7(Sn ( 7 p r 1 ; )i/7 ,(^) i/ ^) fo ry <^ 



From this, we define 



G(p,B,5) := (^J2(B- _2_(5p7-i)i/7). 

Then, f/ r satisfies (12301) and (ET211 if and only if 

G(p P , 5+, 5 + )(l + o; + - w») = m+ (2.24) 

G(p r ,5^,5^)(l + uj- + w*) = (2.25) 

for given by (|2.20p . Solving (|2.24p for w*, and plugging it into (|2.25p . we get 

H(p r ,u) :=G(p r ,B ,S )(2 + uj- ™° ) - = for W = ^++^_. (2.27) 

G(p r ,B£ ,Sfi) 

From (j2"7L3|) and ([2T20]) . we have H(p h 0) = 0. Also, a direct computation using ([H]) and ([2"7T2j) 
yields 

and this implies d p H(pi,0) < 0. By the implicit function theorem, we can choose a constant 
uq > small depending on C/j so that for any u) G [—u>o, u>o], there exists unique p T {uj) satisfying 
H(p r (uj),uj) = 0, and such p r (u) is C 1 with respect to w. Moreover, we may adjust ujq to satisfy 

1T~ = _ ^tIf > for a11 u G [~ w O,^o]- 

du o p H 

Once p r (w) is obtained, then cj* is given by (|2,26p , □ 

Remark 2.2. From Lemma [Kl\ we have that ifui± satisfy (I2.18P with a < ^ for loq in Lemma 
00 tfien 



| U+ -U+\ + \U~ -Uf\ + \u*\< Co (2.28) 
/or a constant C depending only on Ui and 7. 
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2.2. Problems and main theorems. The following is the main problem of this paper. 

Problem 1. Fix oj± satisfying f|2. 18[) and h± satisfying (|2.15p - ([2.17p for a G (0,u;o]. Then, 
find a function go and a vector valued function U = (u,v,p,p) in M to satisfy the following 
properties: 

(i) U is a weak solution to f)2. 1 j> - fj2.3j) in M in sense of f|2.6[) .• 

(ii) U is in C 1 inAf\T D forT D := {(x,g D (x)) : x G R}; 

(iii) U satisfies the R-H conditions (|2.11|) on Tjj; 

(iv) ^ < 1 injV for \u\ = \Ju 2 + v 2 and c defined by (1 1 . 3 [) .' 

(v) U satisfies the slip boundary condition (|2.19p on dJ\f; 

(vi) Setting Af + := M n {y > go{x)} and J\f~ := Af D {y < go(x)}, 

r™oo ^ ~ U fh°°{rt + \{ X >-R}) = ^ W U ~ U r\\L°°(tf-\ {x >-R}) = 0. 2g 

i™o ~ ^ + HL-(A7+\{x<fl}) = ^ W U ~ U ^Wl°°(M-\{x<R}) = 



where U r is as in Lemma \2. 1\ 
(vii) u > and p > ZioZd in N '. 

In order to find a solution to Problem [IJ we use weighted Holder norms. For a connected 
open set f2 C M 2 , let T be a closed portion of the boundary of £1. For x = (x, y),x = (x, y) G £1, 
set 

d x := dist(x, T), d x & := min(d £C , d^). 
For fc£l,a£ (0, 1) and m G Z+, we define 

}2 sup |D^u(aj)|, [u] m ,a,n-= J2 SU P ] _ \* 



\u\\m,o,n 



\u\\ (k ' r) 



0<|/3|<m p|= m 



sup(d a; ) max(l ' 3|+fc ' 0) | J D /3 n( : r)| 



(*,T) ._ V- smaK(m+ a +fc,0) IgMg) ~ D P u{x)\ 

m,a,n-- . S ^P IVJ , _ ~, a 

I II II II I r i ii n(fc,T) n n(fc,T) . r ,(fc,T) 

\u\\m,a,a ■= \m\m,0,n + [%,<.,!i, IMIm,a,n := IML,0,n + Nm.a.n 



where we write = dtd^ 2 for a multi-index /3 = (/5i,/3 2 ) with $ € Z+ and |/3| = fa. + fa. 
Cf^'^JQ) denotes completion of the set of all smooth functions whose || • ||^'J^-norms are finite 

in the norm || • ||~'Jq- Similarly, C m ' a (£l) denotes completion of the set of all smooth functions 
whose || • || miQj n- n orms are finite in the corresponding norm. Hence, Cl^'-^^O) and C m ' a (£l) are 
Banach spaces. For a vector valued function Q = (qj)jLi, let us set 

N N 

wnw - V s II ii nr)ii( fc ' T ) - V ii iK fc ' T ) 

\m\\m,a,n ■— \\Qj\\m,a,n, \\W\\m t a,Cl -~ 2-~i \W\\m,a,Q- 
3=1 3=1 
The following theorem indicates that if a is sufficiently small, then Problem Q] has a solution. 

Theorem 1. Fix a G (0,1), and fix u± and h± with satisfying (|2.18p and (|2.15p - ()2.17p . Then 
there exist constants Cq > and o~q G (0, cjq] depending on Ui, 7 and a so £/taf /or any a G (0, o"o], 
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there exists unique solution U = (u,v,p,p) of Problem\J\ with a contact discontinuity y = go{ x ) 
satisfying the following estimates: 

(i) \\gD\\i,a,R < C a; 

(ii) lim \\g D - co*T]\\i, a M\[-R,R} = for r] defined by (|2.14j); 

(iii) For Af^ as in Problem\^vi) , U satisfies 

, .,(-tx,Th, urn) „ ,,(-a,r h ur D ) „ 

u ~ Uf\\\ cl + + U-Ur ) '.' - ' < Coa; 

(iv) [7 converges to the asymptotic states JJ\ at x = —oo and U r at x = oo in the following 
sense: 

lim \\U-Un ( ' a i h + ^ TD) = lim \\U-U r \\ { ~ a i h -^ D) = 0, 

, ,,(-a,r h , urn) „ „(-a,r h ur n ) 

lim [/ - C/+ ' *+ ^ = lim U - U~ ) ' h " = 0; 

(v) p > -j7j min(p^~, p ; ~) > and u > min(-u^~, uj~) > hold in Af . 

Theorem 2. Fix uj± and h± with satisfying f|2. 18|) and (I2.15p - (|2,17p with a G (0,<7o] for ctq in 
Theorem [TJ Then, there exists a constant C\ > depending only on U\ and 7 so that for any 
(U,go) satisfying all the properties in Theorem^ there hold 

W-UoW^KCta, (2.30) 

\\g D - u*r]\\ L 2(m < Cia (2.31) 



for U (x,y) 



(l-rj(x))(uf,0,pi,pl) + rj(x)(u+,0,p r ,p+) for y > g D (x) 
(1 - 7/(x))(uf ,0, pi, p z ~) + r/(x) (1*7,0, p r ,Pr) for y < g D (x) 



From (12. 11 H . the tangential of To is parallel to the velocity of flow on both sides of r_o. So 
if we use the Euler- Lagrange coordinate transformation, then Tjy becomes a fixed flat boundary 
in the new coordinates while Tp is a free boundary to be found simultaneously with a weak 
solution U of the Euler system to solve Problem [TJ Moreover, due to the conservation of total 
mass flux in x-direction in A/ r± , Fh ± become flat in the new coordinates. 

Let (go, U) be a solution of Problem [lj and for (x, y) G Af, define a transformation T by 

T= (Ti,T 2 ) : (x,y) ^ (x, I pu(x,t) dt - m ). (2.32) 

Jh-(x) 

By (EU, (gJU) , Problem rjjii),(v) and (|2TTD . we have 

-^T 2 (x,g D (x)) = —T 2 (x,h ± (x)) = for any x G M. (2.33) 

By ([233]) and (|iT25j> . we obtain 

T(Af) = Ex (-m^,m£) =: A/", 
T(T D ) = {(X,Y) £ R 2 : 7 = 0} =: T , 

(2.34) 

r(r ft _) = {(x,y) £ 1 : y = -m Q } =: r~, 

T(r h+ ) = G M 2 : Y = m+} =: T+. 

For convenience, we use the notations of 

Af+ =Mn{(x,Y) g m 2 : y > 0}, W~ = A/"n{(x,r) g m 2 : Y < 0}. 
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Since det DT = pu, T is invertible if U satisfies Theorem QJv). Then V = (u, v,p, p)(X, Y) given 
by V = U o T^ 1 is well defined, and it becomes a weak solution of the following system in M: 

dx(—)-9y(-) = (2.35) 
pu u 

d x (u+^)-d Y (—) = (2.36) 
pu u 

d x v + 8 Y p = (2.37) 

2 (U + " } + = j Bfl - for y < (=: B0) - (2 ' 38) 



If a weak solution V of (12T351) — ff2T38|) is in C 1 (7V r± ) n C°(7V ± ), then the integration by parts 
combined with the fact that v = (0, 1) is a unit normal on T° yields the R-H condition 



[^]r° = [p]r° = 0. (2.39) 

Conversely, we can show that if V G C 1 (N' ± ) n C°(A7±) is a weak solution of (12351 

then there exist inverse Euler-Lagrange transformation T so that C/ = V o T -1 satisfies all the 

properties of Problem [TJ 

Lemma 2.3. Fix two C 1,a functions h±. There exists a constant 5\ > depending only on Ui 
and 7 so that ifV = (u,v,p,p) satisfies the following properties: 

W II v - ^ + IIl-(aa+) + \\v - IIl°°(jv-) < 



(ii) Visa wea/c so/utton o/ (I2~35l) - (l2~38]l . and V is in C 1 ^) n C ^); 

(iii) V satisfies the slip boundary condition 

^(X, -m ) = ti_{X) on T, ~(X,m£) = h' + {X) on T + (2.40) 

then, 

(a) the transformation 

f Y 1 

% : (X, Y) ' ^ (X, / —(X, t)dt + /»_(*)) =: (x, y) (2.41) 

is well defined, and satisfies % = T~ x for T defined by (|2.32p ; 

(b) C7 = V o%~ 1 satisfies (i)-(v) of Problem^ with Fry given by 

T D = T- 1 ({Y = 0}) = {(x,y):y= f — (x,t)dt + h^(x)}. 

J-m- PU 

Proof. By using (I2.35p . (I2,40p . (12.41 H . and choosing 5\ small, one can directly check that -DT = 
(DT)~ X . Also, we have T({y = — tUq }) = {y = h-(x)}. This proves (a). 

Using (|2.35p - (12.38p and (I2.40p . one can directly show that U satisfies properties (i),(ii),(iv),(v) 
of Problem [T] if we reduce Si > depending only on Ui and 7. By the R-H condition (I2.39p . U 
satisfies p(u - u)[u]r D = for u = (u,v) where v indicates a unit normal on Tp, and this implies 
either u • v = or [u]r D = on Fd- From the definition of Ui in ()2.13p . we have u^~ 7^ uT ■ we 

further reduce 5\ to satisfy 5\ < — 10 "' ■ . Then we get |[«]r D | > 0. So we have u ■ v = on 
T D . ' □ 

By Lemma |2.3[ Problem [1] can be replaced by the following problem: 
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Problem 2. Fix u)± satisfying (|2.18p and h± satisfying (12.15p - (12,17p for a G (0, cjo] where loq 
is as in Lemma PO . Then, find a vector valued function V = (u,v,p,p) in M satisfying the 
following properties: 



(i) V G C 1 ^) n C ,0 (A^ ± ) is a weak solution of (I2T35I) - (12371) . and satisfies (12351) in N; 

(ii) V is subsonic in M in sense that (u 2 + v 2 ) < c 2 for c defined by (|1.3p .- 

(iii) V satisfies the slip boundary condition (|2.40p .- 

(iv) V converges to the asymptotic states Ui at x = — oo and U r at x = oo in the following 
sense: 

lim \\V - U t \\ Lo0 (^+\i x >-R}) = l im II 17 ~ ^Fllz°°(Af-\{X>-.R}) = °> 

K— »oo H— >oo " 

lir^ II V" - J7 r ||i°°(AA+\{x<ii}) = II V — U~\\ L oo(N-\{ X <R}) =0; 

(v) u > and p > ZioZd in jV. 

Hereafter, we write (X,Y) as (x,y). If V = (u,v,p,p) is in C , then using (|2.35p . (|2.37p and 
(|2.38p . we can rewrite f)2.36j) as 

8J—) = 0. 

So if V is a solution to Problem [21 then we have 

--:S ). (2.42) 



p I Sj" for y > 



[5q for y < 

From (|2.35|) . we expect that there is a function tp satisfying 

v 1 

ip x = -, </? y = — (2.43) 
u pu 

so that we can rewrite (I2.37P and (|2.38p in terms of tp Xl tp y , So and Bo as in [3]. If so, by (|2.40p . 
tp should satisfy 

ip = h^ + k^ on r~, tp = h + + k + on T + (2.44) 

for some constants k±. First, we define tp corresponding to the asymptotic states U\ and U r as 
follows: 

(-*?(=: <pf) fory>0 (^ T +uJ=:<p+) for y > 

|^H^) forn<0 [ + w,(=: ^ ) for y < 

Then, we have 

tpi(x,y) = lim h±{x) and tp r (x,y) = lim h±{x) on T 1 * 1 . 

From this, we choose fc± = in the boundary condition (|2.44j) . 

From (|2.42p . we have p = SoP 7 , and by plugging this into (p 2 ■ (|2.38p ) and using (|2.43[) . we 

get 

&(p,Dtp) = (2.46) 

for © defined by 



^ (i> . q) : = Bop 2 - ^-SofP +l - ^4- for 9 = (<?1> ^2) G M 2 . (2.17) 

7-1 2g£ 
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A direct computation shows ®(pf,Dtpf) = and d p <8(pf ,Dtpf) = -pf{{cf) 2 - {uf) 2 ) < 
in M ± , so we can choose a constant 62 > depending only on Ui and 7 so that for any 
(x,y,q) G (J\f + x B$ 2 {Dtp~^)) U (A/ -- x Bs 2 (DtpJ~)), there exists unique p = p(x,y,q) satisfying 

®(p(x,y,q),q)=0, (2.48) 

and such p(x, y, q) is continuously differentiable with respect to q. We write as p(x, y, q) rather 
than p(q) because Bq in (|2.38p and So hi (12.42|) are piecewise constant functions. 

Let (p G C°(Af) fl C 1 (Af ± ) be a function satisfying \\<p — ( ^z + IIc 1 (a7+) + 11^ — ¥ , rllc 1 (j\F') — 
then there exists unique y, Dip) satisfying the equation 

<d(p(x,y,Dip),Dip) =0 in M ± . (2.49) 

For such p(x, y, Dip), we use (|2.42p and (12.43P to express u, v and p as 

u = , v = -^-, p = Sop 1 (2.50) 

pp y ptpy 

so that the equation (|2,37p can be rewritten as 

div(A(x,y,D<p)) = in N (2.51) 

where A(x,y,q) = (A\, A2)(x, y, q) is given by 

A l (x,y,q) = — — 1 . , A 2 (x,y,q) = S (x,y)p' y (x,y, q) (2.52) 

with 

W A j\\c*mxB 3S2/ t{D V +)) + W A ^C^-XB^ /A {D V -)) ~ Ck ( 2 - 53 ) 

for each k G N and j = 1, 2, where the constant depends only on Ui, 7 and fc. 

Next, we consider the R-H condition for 99 on T°. Rewriting (|2.39p in terms of tp, we get 
[Vx]r = [A2(x,y, D(p)] r o = 0. In particular, we may rewrite [<^a;]r = as [np]-po = k* for a 
constant fc*. Furthermore, we choose k* = by continuity of </?j and (p r across T°. Because, we 
will seek a solution ip of (|2.5ip so that ip converges to cpi at x = —00 and to tp r at x = 00. So 
we get 

[<p]TO = [A 2 (x,y,D<p)] T o=0. (2.54) 
Now we consider the boundary value problem (|2.5ip . (I2.44p with k± = and (I2.54p . 

Theorem 3. Fix a G (0,1), and fix ui± satisfying (|2.18p and h± satisfying (|2.15p - (|2.17p . Then, 
there are constants C2 and o~\ G (0, ljq] depending on Ui,j and a so that wherever a G (0, <n], 
the boundary value problem (|2.5ip and (|2.44p with k± = /ias unique weak solution ip G H} (J\T) 
satisfying the following properties: 

(i) 99 is in C°(A7) n C^ATi) n C 2 ^), hence <p satisfies (l2~5Tj) in AA ± , and (12341 on T° 
pointoise/y; 

(ii) The equation (|2.5ip is uniformly elliptic in N ; 

(iii) (/? satisfies the estimate 

n + n(-i-a,r + ur°) n _||(-i-a,r"ur°) 

(iv) ip converges to <pi at x = —00 and to (p r at x = 00 in i/ie following sense: 

r 11 +n(-i-aT + ur°) r n _n(-i-a,r-ur ) n 
Jim ||p - <p] \\ 2 ^ + ^_ R} = Jim ||y - ^ II W"\{*>-K} = 0. 

r n +n(-i-a,r+ur ) ,. n _||(-i-a,r-ur°) n 
J™ \\<p - y r \\ 2:QM+ \ {X < R} = hrn^ \\p - <p r ||^ aiAr - Ux <4 = 0; 

(v) dyip > — holds for some m* > in N where m* > depends only on Ui and 7. 



(2.56) 
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We first prove Theorem El then prove Theorem [T] and Theorem [2j 



3. Proof of Theorem [3] 

In order to prove Theorem^ we need to prove unique existence of a weak solution <p £ H^ oc {J\f) 
to the boundary value problem 

div(A(x,y,Dip)) = in TV, (3.1) 
ip = h— on r~, tp = h + on T + (3.2) 

for A defined by (|2.52p . Furthermore, the weak solution (p is required to satisfy the additional R- 
H condition ([2341 on T°. If the equation (J33J) is strictly elliptic and <p £ C 1 (77±)nC 2 (A/' ± ), then 
the weak Harnack inequality implies that <p is continuous across T°, and [A2(x, y, Dip)]-po = 
easily follows from the integration by parts because the vector (0, 1) is unit normal of T°. 
Therefore the key point to prove Theorem [3] is to show that the boundary value problem (|3.ip . 
(|3.2|) has a piecewise C 1 weak solution. For that purpose, we employ results from |12j . 

We prove Theorem [3] in two steps. First, we formulate a linearized boundary value problem 
where coefficients of an elliptic equation in the boundary problem are piecewise C a , and apply 
the result of [12] to weak solutions of the boundary value problem. Main difficulty in the first 
step would be uniform L°° estimate of weak solutions in unbounded domain Af because of two 
different asymptotic states ipi and <p r at x = ±oo. Then, we use a fixed point theorem to prove 
Theorem [3l 



3.1. Linearized boundary value problem. Define a smooth connection from ipi to ip r as 
follows: For r/ defined by (|2.14p . we set 

<Po{x,y) = (1 - r]{x))pi{x,y) + r](x)ip r (x,y) (3.3) 

for <pi and ip r defined by (|2.45p . By (|2.28|) . if we choose a small, then A(x, y, Dpo) is well defined 
by (|2.53p . Since div(A(x,y, Dpi)) = 0, (|3.ip is equivalent to div(A(x,y, Dip) — A(x,y, Dtpo)) = 
divFo in M for 

F o :=A(x,y,Dipi)-A(x,y,D<p ). (3.4) 
To ensure that Fq is well defined, we let 

a < min{o;o, ^} =: cr fl (3.5) 

for C in (^28]) . From §ZM§ , <l23U|) . (l232|) . (I3T3]) and the definition of F in ([33]), we easily get 

„ . N [(0,0) for x < -1 , . 

F (x,y)= ' ; " (3.6) 

I {u,pi — p r ) tor X > 1. 

Moreover, ([2T28T) implies that F £ C°°(7?+) n C°°(JF) satisfies the estimate 

ll^ollc-fc^) + \\Fo\\ c >'Q7=) - Ck<T ( 3 - 7 ) 
for a constant depending only on J7;,7 and k for each k £ Z+. For a fixed function </>, set 

I' 1 

af{x,y):= D qj Mx,y,Dip + tD(cP-ip ))dt (3.8) 

J 
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for Ai(x, y, q) defined by <|2.52 j> . ip solves the boundary value problem of (13. ip and (j3.2|) if and 
only if ip := ip — <po solves 
2 

y dMfdjiP) = divFo in M, 

;j: 1 (3-9) 
ip = h— + 1 + (1 — rf)U- =: g- on T~ , ip = h + — 1 — (1 — rj)u)+ =: g+ on F + 

where d\ and 82 denote J| and -g- respectively. 
Fix a G (0, 1), and define an iteration set /Cm by 

/Cm := W e C7^(A^)nC°(AT) : ||^-^lfc 1 ^^ ) + llv-vrll£? r " UI<,) < *M (3-10) 

for constants M > 1 and a G (0, to be determined later with Ma < for 82 in (|2,53p . 
Then, Km is a convex and compact subset of the Banach space C^_"( 2 (A/'+) n Cv^s 2 (AA~ ) (see 
Section 5.1 of [1] for details). For a fixed G /Cm, consider the following linear boundary value 
problem: 

2 

^ ^(ag^) = divFo in jV, (3.11) 

ip = h- + 1 + (1 - rj)uj- =: g- on r~, ^ = h+ - 1 - (1 - 7/)w + =: on T+. (3.12) 
The following lemma is essential to prove Theorem [3l 

Proposition 3.1. Let 82 be as in (|2.53p . Fix a G (0,1). Then there exist constants C > 
and eo G (0, ^p] depending only on f/j,7 and a so i/iai if Ma < £0, then for any (ft G /Cm, 
boundary value problem of (|3.1ip and (|3.12p /ias unique weak solution ip G Hl oc {N) H C (A/") 
mi/i satisfying 

w&r^ + ii*ii&°- r_ur "* (3.i3) 

< C([|- F o[|l,a^V+ + ll-fblll,^- + l|ffl|l,a,AA + H-Dfllka^)' 
ana" Jim [Hlg^J^ = Jim = (3.14) 

where g is defined by (j3. 15j) . 

In order to prove Proposition 13.11 we set a linear boundary value problem in a bounded 
domain as follows: Let xo be a function satisfying 

Xo(y) = l 1 fOTy -:3 L ', 0< Xo (y)<l for y6l4.fl- 
[0 for y > z z 

xo(y) < for a11 y e M and llxollc 2 (K) < c("io) 

for a constant C(m^) depending only on for in (|2.2U|) . and define a function 5 by 

g(%,y) ■= g-(x)xo(y) + g+(x)(i - xo(y)) (3.15) 

for given by (|3.9p . Then 5 satisfies the estimate 

IMIl,a,JV ^ C(ll<7-||l,a,R + ||5+||i )Q! ,r) 

for a constant C depending on Xl\ and a. For a fixed constant R > 10, set Mr := M H {(x, y) G 
M 2 : \x\ < R + 2}, and let AAr be a convex and connected domain satisfying Nr C A//?, C Nr+i 
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where 8Mr is a simple closed smooth curve. We consider the following boundary value problem 
in a bounded domain Mr: 

2 

diiafdj^) = divF in Mr, (3.16) 

iP = g on OMr. (3.17) 
(|3.16p and (|3.17p has unique weak solution i/jr in ^(Mr). We claim that ^ = lim ipR is unique 

R— ¥OD 

weak solution of (l3~TTT) and (pTL2l) with satisfying (f3TT3l) and ([37HD . 

Lemma 3.2. Lei 52 &e as in (|2.53p . There exist positive constants C > 0, £i G (0, ^l 2 -] and A 
depending on XJ\ and 7 uraia C depending on a in addition so that if Ada < e± then, for any 
<f) G Km , the coefficient matrix (a^J defined by ()3.8p satisfies the following properties: 



W A|^| 2 < E^i^foy)^ < ll^l 2 »n for any £ = € 
(ii) (a^) 2 =1 is symmetric, that is, a^(x,y) = a^(x,y) in M; 

m - 4? ) ii;;:^ ur * ) + ii^' - "Sr , ii^- ur0) * cm, /or M = 1, 2. 

Proof. Set c 2 = 7 So/) 7 " 1 . Then a direct computation using (I2.47jl and (|2.48j) yields = grrj^p 
and d q2 p = p "i" + X - for |u| 2 = -a 2 + 1> 2 , and this implies 



a 4 c 2 -n 2 /? 2 c 2 n|it| 2 /?c 2 m; 

fl^Ai = " c 2_i M |2 ' 8 >n A * = c 2 _ up 1 Vi = V2 = ~ ^rn p - ( 3 - 18 ) 



.2„4„2 



This proves (ii). Using and (l3~T8]l . we can easily show that ag"ag" - (a£ l) ) = %f^k > 

f 1 \ f(pt> u t>Pl) for y > 

for (fiuuupi) = <._'_ , , ,sowe have 

[(.Pi , u i ,pi) fo r y < ° 

A |€| 2 < X>SfW)^< fl£| 2 

for all (x, y) E M and £ = (^1,^2) G M 2 where Ao > is a constant depending on U\ and 7. Since 
A(x,y,q) is smooth with respect to q near D(pi, one can choose e\ sufficiently small depending 
on U[ and 7 so that if Ma < e% in the definition of ICm, then we obtain (i) of Lemma 13.21 for 
A = jg 1 . (hi) can be easily checked from (|3.8p and f|2.53|) . □ 

Remark 3.3. Lemma \S \2( ii) is a necessary condition to apply the result of |12j . See Theorem 
1.1 of [12J for details. 

Proposition 3.4. There exists a constant C depending only on Ui,j and a so that if Ma < E\ 
for £\ in Lemma \3.2\ and R > 10, then for any <j> G ICm, the boundary value problem of f|3. 161) 
and (|3.17p has unique weak solution ip G H 1 (Mr) n C°(Mr) satisfying the estimate 

(-1-a.r+ur ) I, ,|,(-i-Q,r-ur ) 

< C( 11-^0111,0:,^+ + H-fb||l,a,AT- + llfl'l|l,a,JV + \\Dg\\ L 2 {Af) ) 

where we set M^ := {(x, y) G Mr :y>0}, M R := {(x, y) G Mr :y<0}. 
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Proof. (Step 1) Since the equation (13.16P is uniformly elliptic, the boundary value problem of 
(I3.16P and fj3. 17f) has unique weak solution ip in H (Mr). Also, the Harnack inequality (see [10] 
Theorem 4.17, Corollary 4.18]) implies that ip is continuous across T° C\Mr. 

By Proposition 3.2 of [12], we can choose ri(< R/W) small and C depending only on t/j,7 
and a such that for any Zq G A/j? D {(x,y) : \x\ < R — 2r\, \y\ < y^}, there is a continuous 
piecewise linear function LS Z ^ satisfying 

|^(x) - L (Zo) (x)| < C7/i |x - Z \ 1+a for all x G B ri (Z ), (3.20) 

L (^») G ^1(^7+) n c^AF), 11^^)11^^) < C/xo (3.21) 

with 

MO = IMIl^A/'r) + ll- F o||l,a,Ar + + 11-^0 1| l,o,Af- + IIS'l|l,aX- (3.22) 

Then, by using the method of the proof for Theorem 1.1 in |12| . we obtain 

HV'lli > a^ /B n{|y|<a} + II^Hi^ 8 n{|y|<a} ^ (3 ' 23) 

for a constant C depending only on a. Differently from Theorem 1.1 of [12], we note that the 
regularity of ip is not weaker than C > a because T° is flat thus in C°°. By Theorem 8.33 in [9] 
and LemmaE^l we also have H^L ^+ n( , , + ^- n r|,,i^nT. < Ca*o for a constant 

' > JV 4fl/5 l llyl^20- r ' > JV 4ii/5 'U»l ^ 20 J 

C depending only on t/j,7 and a, then combining this estimate with (|3.23p yields 

11^1,^ + W;^,^^ ( 3 - 24 ) 

For a fixed point Xo G M^ 2 , let 2d* := dist(Ao,r + U T ). It suffices to consider the case of 
24 < n. Set WWfa) := ^o+d^-^Xo) for ^ £ B (Xo) . = {rj £ Bi(Q) . Xq + ^ £ A^J. 
Then the standard elliptic interior estimates yield 

ll^ (Xo) ll 2aB C*o) < C(U\\ la ^. , + ll^o||i >a ^+) ; (3-25) 

and by scaling back and combining (|3.25p with (|3.24p . we get 1 UF ' < C/^q- Repeating 

2,a,M R/2 

the same argument for points in A/"^ 2 j we S e ^ 

I^H^ r+Ur0) + Wtl-ZT^ < ^o- (3.26) 

(Step 2) In order to finish the proof, it remains to estimate ||V'IIl° o (a/'h)- Set u := ip — g then 
u G Hq(Mr) satisfies 

/ afdjudiC = [ F ■ D( - afdjgdiC for any ( 6 ^o(A^). (3.27) 

By (|3.6p . we have 



F -D(= Fo . DC+ ( pi - Pr ) dyC= Fq . DC 

Mr JAf R n{\x\<2} JAf R n{x>2} JAf a n{\x\<2} 

Plug C = u into (|3.27p . Then, by Lemma l3.2f i) and Holder inequality, we get 

/ \Du\ 2 < C{\\F4 2 Loo[M) + \\Dg\\ 2 L2(M) ) =: Ml (3.28) 

where C depends on Ui and 7 by Lemma 13.21 but independent of R. Fix a positive constant Iq 
for to < < 1- For each zo G T D Mr, set Qi (zq) := Mr fl {|£ — xq| < Zo}- Since u = on 
<9A/r, the Poincare inequality provides IMI l 2 (q ;o (z )) — C\i\ for all zq G A/j where C depends 
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on Iq but independent of R. Then, by the method of Moser iteration, we can find a constant C 
depending only on Ui and 7 so that u satisfies 

\\u\\l^(M r ) < Cfix. (3.29) 
Combining (^29]) with (I336|) . we finally obtain ([535]) . □ 

3.2. Proof of Proposition [3TT1 Now we can prove Proposition 13.11 easily. 

We choose 8q = E\ for E\ in Lemma 13.21 For each m G N, let ^( m ) be unique weak solution 
of (13.161) and (13.171) in M m -\-20 with satisfying (13.19j) . Then we can extract a subsequence, still 
written as {ip( m '}meNi so that the subsequence converges to a function tp* G Hl oc (N) in the 
following sense: for any R > 

(i) i/;( m ) uniformly converges to tp* in ATr; 

(ii) converges to ^* in C 1 in A/"^" and Nr', 

(iii) ■i/ ;( ' m ' ) converges to tp* in C 2 in if+ and K~ for any if + CC Af + , K~ CC A/" - . 
Also, ip* satisfies the estimate 

H / * || (-i-a,r+urO) || ^..(-l-a.r-ur ) 

II ^ H2,a,Ar+ H2,a,AA- (3 30) 

< C(H-fb||l,a,Ar + + ll-folll.a.AT- + llslll.aA/" + II D 9 Hi 2 (AT) ) 

for C same as in (|3.19p . We claim that tp* is the unique weak solution of (13.111) , (|3.12p satisfying 
the estimates (|3.13|) . (|3.14p in Proposition 13. H 

For each constant R > 10, let xr be a smooth function satisfying 



XR(x,y) = { . - l , \D X r\<W. (3.31) 
10 if \x\ > R — ~, 

Since each is a weak solution of (|3.16p and (|3.17p . by the dominated convergence theorem 
and (^271) . u* := tp* - g £ H} oc (N R ) satisfies 

afd j u*d i (u*x 2 R)= [ F -D(u* XR )-afd 3 gd t (u* XR ) (3.32) 
■ J JN 

from which we get 

f \DtP*\ 2 dX <Cm (3.33) 

for hi in (|3.28p by (|3.29p where the constant C depends on U\ and 7 but independent of R. 
Since R can be arbitrarily large, we get 



—^00 



lim / \DtP*\ 2 dx = 0. (3.34) 

/Mn{\x-Q\<2} 



Since g± G H 1 ^), we have Jim Jj a ._Q|< 2 m^)\ 2 dx = Jim ij a; _Q|<2 Iff-O^i m o )\ 2dx 

g{x, m+) - J™° ^*(x, t)cft for y > 



0. Then, expressing ip*(x,y) as tp*(x,y) 
easily show that 



, we can 

g(x, -m ) + J_ m - ip*(x, t)dt for y < 



and this implies that 



lim / (V>*) 2 dX = 0, (3.35) 

IQI^oo J\ x -Q\<2 

IQ^OO H^*Hi°°(|x-(5|<|) = - ( 3 - 36 ) 
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by (13, 6h and [HI Theorem 8.17 and 8.25]. Repeating the argument of (Step 1) in the proof of 
Proposition 13.41 with using (12.15P and (13.6p . we can show that ip* satisfies (13.14p . 

Suppose that Vr and ip2 are weak solutions to (I3.11j) and (I3.12j) . and that they satisfy (13.13P 
and (|3.14p . Then, u := rpx — ip2 satisfies Jj^a^djud^uxfydX = for any R > 10, and from 
this and (|3.14p . we get 

/ \Du\ 2 < C lim / \u\ 2 = 0. 

J AT R ^°° JAfn{R-l<\x\<R-±} 

This implies Du = in AA ± , therefore we get u = k in Af for some constants k. Since ipi = V>2 
on T^, we conclude that k = thus tpi = ip2 in N . The proof of Proposition 13.11 is complete. □ 

3.3. Proof of Theorem [3J Finally, we prove Theorem [3l 

Fix u>± and h± with satisfying (]2TT5]) - ([2TT71) and ([27TS]) . By Proposition EH if Mcr < ej for 
ei in Lemma |3.2|, then for any <j> G /Cm> the linear boundary value problem of (|3.1ip and (|3.12D 
associated with cf> has unique weak solution ip^> G Hl oc (N) satisfying the estimates f|3. 13|) and 
(|3.14p . We define a mapping X by 

1:^ + tp 

for 990 given by J33J. By (l2~15]l . (l2~T6]l . IBlSjl . (CO]> . (13^1 . (l3T2]l and (13T3L we have 

ii-T-/ |,(-i-a,r+ur°) in-/ ,\ i,(-i-a,r-ur Q ) . ^ 

for a constant C 1, depending only on {//,7 and q. We choose M and <ti by 

M = max{8C*,2}, ^ = min{^, °— , 1} (3.37) 

for cr" in (]3.5h and £i in Lemma 13.21 For such choices of M and o\ , the iteration mapping X 
maps /Cm into itself wherever a < a\ for cr from (|2.15p and (I2.16p . We point out that the choices 
of M and o\ in (|3.37p depend only on t/j,7 and a. 

We claim that X : K>m — > K-M is continuous in C'^ Q j{ 2 (A/* + ) n C*_^ 2 (A/" - ) (see Section 5.1 in 
[4] for the definition of C^°^ 2 ). Suppose that a sequence {4>k} in /Cm converges to <f> G Km in 

CM 2 (AT+) n cM 2 (AF). Let us set ip k := X(<f> k ) - ip for each k G N and ^ := X(0) - By 
(|3.13p . any subsequence of {tpk} has its own subsequence that converges to a function ip* with 
V>* + tpo G /Cm in C^^AF) n CJ^^AF), and such V* is a weak solution to 

2 

^ di^fdjip*) = divF in A/", V* = 5 on r+UF. 
By repeating the argument in the proof of Proposition l3.lt we can show that ip* = vp in N . This 

in C 1 { f/ 2 (N+) n CW 2 (N=). 



implies that tpk converges to ift in C}f,( 2 (N + )nC}^,( 2 (N ). Hence, X : /Cm — > /Cm is continuous 



As pointed out earlier, /Cm is a convex and compact subset of C^^ 2 (A/* + ) n C l ^ 2 {M ). 
Then, by the Schauder fixed point theorem, we conclude that for any given h± satisfying (|2.15p - 
(|2.17p . X has a fixed point <p* in /Cm- By Lemma 13.21 ^ satisfies (ii) of Theorem [3j Also by 
Proposition 13. II and (|3.3p . tp* is a weak solution of (|2.5ip in J\f, and satisfies the equation (|2,5ip 
pointwisely in Af^. Then, for any £ G Cq°(J\T), we have 

I [A 2 (x,y,D<p% ti= f ' A(x,y,D<pl)-DZ+ [ S t divA(x, y, Diffi) = 0, 

Jr° J JN+uAf- 
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so (fr satisfies the R-H condition (12.54p . We may reduce o\ in (I3.37j) further so that (12.55P 
implies (v) of Theorem [3l Then satisfies all the properties stated in Theorem [3j 

Given h = (/i + ,/i_), let <p>^ and ip^ 2 ' be weak solutions of the boundary value problem of 
(|2.5ip and (I2.44p with k± = where cp^ and (p^ satisfy all the properties stated in Theorem 
[3] as well. Let us set ip^) ■= (pW — ip for j = 1, 2. Then ip^ 1 ' — ip 1 ^ satisfies 



(3.38) 



diia^dj^ - ^) + K ( J (1)) - af^d^) =0 in TV, 

- ?/;( 2 ) = on dM. 
By (12321) and (|3T8]), we can rewrite (EOBD as 

E^((a;f )) +6 u -H-w (1) -^ (2) )) = o 

with || &ij || < CMa. We again reduce <7i > depending on Ui and 7 to have 

— |r/| 2 < Yl + - yl^l 2 in -A/" for all 77 = (771,772) 6 M 2 

for some constant A' > 0. Then, repeating the argument in the proof of Proposition 13.11 we 
conclude that = ip^ in M, thus tp( 1 ' = (p^ in M. Finally, we choose C2 = 2M, then the 
proof of Theorem [3] is complete. □ 

4. Proof of Theorem Q] and Theorem [2] 

4.1. Proof of Theorem[T]. First, we choose do = <J\ for <7i in Theorem[3l and fix d± satisfying 
(|2.18p and h± satisfying (|2.15p - (|2.17p for a G (0, o"o]- Let tp be the corresponding solution 
satisfying all the properties stated in Theorem [3j Let (u,v,p) be given by (12.50h from <p with 
p determined by (|2.49|) . then V = (u,v,p,p) satisfies all the properties stated in Problem [2] as 
well as the estimates 

|m 11 v - un[;:^IU = w - "ntT^U = °> m 

lim \\V - r/+||(- a ' r+ur °) - l im \\ V _ tt- ur°) _ 

By (|4. 1 j) . the inverse Euler-Lagrange transformation T given by (|2.4ip is continuously differen- 
tiable in and invertible with the estimates ll^llca,^;^!) < C and a /^±j < C- From 

this, one can directly check that ^jd^(X) < d^(X) < Cd^(X) for any X G AA + where we 
set d£{X) := dist(X,r+ U T°) and d^(X) := dist(T(A"), T(r+ U T )), that is, d\ is a distance 
function in the Lagrangian coordinates, and is a distance function in the Eulerian coordi- 
nates. Similarly, if we set d L (X) := dist(X, F~ U r°) and d E (X) := dist{%{X),%(T~ U T )) for 
X G W~, then we also have ^d L (X) < d E {X) < Cd L (X). 
If we set 

r D :=T(r°) and U:=Vo t 5T 1 , 

then U satisfies all the properties of Problem [1] as well as (iii) and (iv) of Theorem [1] where we 
choose Cq as CC\ for some constant C depending only on Ui,j and a. By the definition of T in 
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(|2.4ip . the contact discontinuity is given by 

f° 1 

T D = {(x,g D (x)) : g D (x) = \ — - — -rdt + h-(x),x G M} 



pu(x, t) 

r m t l 

or = {(x,g D (x)) : #d(x) = h+(x) - / — - — — dt, x E R}. 



(4.2) 

dt t a Ml 

lo pu(x,t) 

By (|2.43p . (|2.44p with fc± = and (|2.45[) . we can express g D as 

/*° 

g D (x) = I _(py(x,t)dt + h-(x) = tp(x,0) = (ip — tpi)(x,0). (4.3) 

Then, (i) and (ii) of Theorem [T] follow from (|2.45p . (|2.55p and (I2.56p . Theorem [T^v) follows from 
Theorem ffly) and (jXMD - 

Now it remains to verify the uniqueness in Theorem [TJ For fixed u± and h±, let U^> and 
be two solutions of Problem [1] with satisfying all the properties stated in Theorem [JJ and 
let T® for j = 1,2 be defined from by (l232|) . and set := o (T^)" 1 in A/". Then, 
each = (u^\ v w\p^\ p^) satisfies (I4.ip . For j = 1,2, if we can find <p^' satisfying 

(J) 1 

'^■^-U'^oSo* inA ^ (44> 

pC*') = fr ± on T ± (4.5) 

then </7^^ satisfies all the properties in Theorem [3] by reducing o"o > if necessary. Then, 
following the proof of Theorem [31 we can easily show that ip^ 1 ' = (p^ in M so = in 
J\f. The proof is complete. □ 

4.2. Proof of Theorem [2j For u± and h± fixed, let {U,go) be the corresponding solution 
satisfying all the properties in Theorem [TJ and let <p* be the solution corresponding to (U^go) 
so that ip* satisfies all the properties in Theorem Then, by (|2.43|) and (|4.2|) . the contact 
discontinuity function go is given by 



g D (x) = h + (x) - (p y (x,t)dt. 
J o 

Then, by (|230|) . (1X151) and (ETT2"]) . we have s£>(x) - w*t/(x) = g+(x) - J Q m ° <9 y (v? - (f )(x,t)dt 
for in Lemma 12.11 and r] defined by (|2.14p , and this implies 

Wqd ~ ^*^||l2 (m) < 2(||0 + || L 2 (M) + m$\\D((p - <A))IIl 2 (a/'+))- 

Then, by (T2TT5D . (T2T26D . q> and (13^31 . we obtain 

(1230]) can be similarly proved by using (12381) . (g32D , (|333l) and Theorem [If v) . □ 

Acknowledgments. The author thanks Gui-Qiang Chen at University of Oxford for moti- 
vating to work on this problem, and Mikhail Feldman at University of Wisconsin-Madison for 
helpful discussion. 
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